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ABSTRACT 


Recent  observational  and  modeling  studies  have  suggested  that  transient 
tropical  perturbations  may  induce  significant  intraseasonal  teleconnections 
between  the  tropics  and  midlatitudes.  We  have  investigated  the  mechanisms  for 
such  teleconnections  using  a  nonlinear  global  shallow  water  model  with  a 
realistically  wavy  basic  state.  The  model  perturbations  were  designed  to  simulate 
tropical  cyclones  with  prescribed  growlh,  decay,  and  propagation. 

The  model  responses  to  the  tropical  perturbations  showed  distinct 
midlatitude  wavetrain  responses  to  the  perturbation.  The  typical  response  became 
large  within  a  few  days  and  remained  as  a  coherent  pattern  for  two  to  three  weeks 
after  the  demise  of  the  tropical  perturbation.  The  response  was  particularly  strong 
if  the  perturbation  propagated  close  to  an  extratropical  jet.  This  propagation 
allowed  Rossby  wave  induction  as  divergent  outflow  from  the  perturbation  crossed 
the  jet's  region  of  high  ambient  vorticity  gradient.  This  initial  wave  response  was 
then  guided  by  the  jet  and  amplified  in  barotropically  unstable  regions  on  the  jet 
flanks.  This  produced  a  persistent  response  downstream  of  the  jet  exit  (e.g.,  over 
the  northeast  Pacific  and  North  America). 
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The  response  showed  a  marked  variation  with  the  season  and  with  the 
tropical  ocean  basin  in  which  the  perturbation  occurred,  apparently  because  of 
temporal  and  spatial  changes  in  the  basic  state  conditions.  The  responses  closely 
resemble  observed  anomaly  patterns  associated  with  transient  tropical 
perturbations,  and  are  dynamically  consistent  with  the  responses  to  more  persistent 
tropical  perturbations  (e.g.,  ENSO  events). 
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I.  INTRODUCTION 


A.  BACKGROUND 

Understanding  climate  has  been  an  object  of  study  throughout  history  due  to 
climate's  large  impacts  on  human  endeavors.  Early  observers  realized  the  importance  of 
the  net  effects  of  the  atmosphere  on  a  continental  or  even  a  global  scale.  Drought,  floods, 
and  famine  have  occurred  throughout  history  in  an  impredictable  manner.  Yet  many 
aspects  of  climate  (the  large  scale,  low  frequency  state  of  the  atmosphere)  showed  distinct 
spatial  and  temporal  patterns.  These  patterns  have  long  motivated  researchers  to  explore 
the  underlying  processes.  These  climate  studies  continue  today,  with  the  objective  of 
making  long  range  atmospheric  predictions  that  would  provide  adequate  warning  of 
adverse  environmental  conditions. 

B.  GLOBAL  SCALE  ATMOSPHERIC  INTERACTIONS 

Walker  (1924)  and  Walker  and  Bliss  (1932)  conducted  seminal  climate  research 
on  global  patterns  of  atmospheric  pressure,  temperature  and  precipitation  anomalies 
associated  with  monsoonal  circulations  in  the  Indian  Ocean  and  the  western  Pacific. 
Walker  and  Bliss  (1932)  identified  a  major  periodic  atmospheric  mass  transfer  between 
the  eastern  and  western  tropical  Pacific  which  has  come  to  be  known  as  the  Southern 
Oscillation.  Walker  and  Bliss  (1932)  detailed  long  range  statistical  correlations  between 
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the  Southern  Oscillation  and: 


•  Hawaiian  Island  precipitation, 

•  western  Canadian  surface  temperatures, 

•  sea  level  pressure  over  the  southeastern  United  States. 

By  correlating  periodic  variations  in  various  atmospheric  fields  at  distant 
locations.  Walker  and  Bliss  (1932)  produced  the  first  study  of  atmospheric 
teleconnections  (Horel  and  Wallace  1981,  Wallace  and  Gutzler  1981).  Bjerknes  (1966, 
1969,  1971)  furthered  teleconnection  research  by  suggesting  that  large  scale  circulation 
anomalies  were  driven  by  anomalous  tropical  air-sea  interactions  that  produced 
anomalous  energy  and  mass  transports  between  the  tropics  and  the  mid-latitudes. 

Bjerknes  (1969)  furthered  Sir  Gilbert  Walker's  Southern  Oscillation  research  by  linking 
the  oscillation  of  sea  level  pressure  anomalies  to  low  frequency  fluctuations  in  sea  surface 
temperature. 

C.  INTERANNUAL  TELECONNECTIONS 

The  bulk  of  atmospheric  teleconnection  research  has  been  focused  on  the 
interannual  teleconnections  that  arise  due  to  the  phenomena  called  El  Nino  and  its 
symbiotic  relationship  with  the  Southern  Oscillation.  Bjerknes  (1966)  first  used  the  term 
El  Nino  to  describe  a  large  scale  anomalous  warming  of  the  eastern  equatorial  Pacific 
from  the  South  American  west  coast  to  the  international  dateline.  During  El  Nino, 
coastal  upwelling  along  the  western  coast  of  South  America  weakens,  easterly  trade 
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winds  weaken,  and  the  sea  surface  temperature  (SST)  gradient  diminishes  between  the 
western  and  eastern  regions  of  the  equatorial  Pacific.  El  Nino  events  generally  occur 
every  2-7  years  (Rasmusson  and  Carpenter  1982).  Contrasting  events,  sometimes  called 
La  Nina  events  also  occur.  In  these  events,  the  anomalies  seen  during  El  Nino  events 
occur  in  reverse  (Bradley  et  al.  1987). 

Bjerknes  (1 966,1 969,1 972)  was  the  first  to  recognize  the  link  between  El  Nino 
events  and  large  scale  extratropical  atmospheric  circulation  anomalies.  He  related  these 
circulation  anomalies  to  changes  in  the  strength  of  the  Hadley  and  Walker  circulations. 
Horel  and  Wallace  (1981)  described  the  teleconnections  associated  with  El  Nino,  La 
Nina,  and  the  Southern  Oscillation,  emphasizing  the  vertical  structure  of  the  atmosphere 
and  interactions  between  the  tropics  and  the  middle  latitudes.  They  studied  the  Northern 
Hemisphere  winters  between  1951  and  1971  during  which  several  El  Nino  and  La  Nina 
events  occurred  and  obtained  significant  correlations  between  the  interannual  variability 
of  mean  tropospheric  temperature,  200  mb  heights,  SST,  and  many  other  variables.  In 
particular,  the  correlations  of  tropical  Pacific  SST  with  Northern  Hemisphere 
geopotential  heights  exhibited  well  defined  telecoimections  in  the  Northern  Hemisphere. 
The  anomaly  patterns  reported  by  Horel  and  Wallace  (1981)  are  consistent  with  the 
northern  winter  teleconnection  patterns  discussed  by  Walker  and  Bliss  (1932)  and 
Bjerknes  (1966,  1969). 
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Wallace  and  Gutzler  (1981)  also  studied  Northern  Hemisphere  winter 
telecoimections  using  sea  level  pressure  and  500  mb  geopotential  height  analyses  for  a  15 
year  period  and  found  distinct  telecoimection  patterns  not  only  in  the  Pacific  but  in  the 
Atlantic  as  well.  They  described: 

•  a  North  Atlantic  oscillation, 

•  a  North  Pacific  oscillation, 

•  a  Pacific/North  American  (PNA)  pattern. 

The  PNA  pattern  is  an  arc-like  arrangement  of  four  tropospheric  circulation  anomaly 
cells:  an  anticyclonic  cell  near  Hawaii,  a  cyclonic  cell  over  the  North  Pacific,  an 
anticyclonic  cell  over  western  North  America,  and  a  cyclonic  cell  over  the  western  North 
Atlantic.  This  pattern  is  very  similar  to  one  found  by  Horel  and  Wallace  (1981)  (Figure 
1).  Wallace  and  Gutzler  (1981),  Hoskins  and  Karoly  (1981),  and  others  suggested  that 
these  teleconnection  patterns  represent  standing  oscillations  in  the  tropospheric 
geopotential  height  field.  These  oscillations  strongly  resemble  the  standing  wave  pattern 
foimd  in  many  atmospheric  models  forced  by  tropical  disturbances  (e.g.,  Hoskins  and 
Karoly  1981). 

D.  INTRASEASONAL  TELECONNECTIONS 

Teleconnections  may  occur  on  intraseasonal  time  scales,  apparently  due  to 
anomalous  forcing  occurring  on  intraseasonal  scales.  Nitta  (1987)  correlated  tropical 
western  Pacific  SST  anomalies  and  convection  anomalies  with  anomalies  in  mid-latitude 
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tropospheric  winds.  He  found  that  the  intraseasonal  response  to  anomalously  warm  sea 
surface  temperatures  was  more  than  twice  the  interannual  response.  Chang  and  Lum 
(1985)  studied  the  relationships  between  the  northern  hemisphere  winter  east  Asian  - 
North  Pacific  jet  and  nearby  tropical  convection.  They  found  a  strong  positive  correlation 
between  convection  anomalies  and  jet  strength  anomalies  with  no  apparent  time  lag. 
However,  they  suggested  that  the  tropical  convection  forced  the  strengthening  of  the  jet 
via  a  process  in  which  increased  upper  tropospheric  divergent  flow  from  the  convective 
region  experienced  Coriolis  turning  and  intensified  the  jet  flow. 

This  extratropical  jet  response  to  increased  tropical  convection  has  also  been  seen  in  the 
southern  hemisphere  (Hurrell  and  Vincent  1991). 

Harr  and  Elsberry  (1995)  conducted  an  observational  study  which  showed  that 
western  Pacific  tropical  cyclones,  especially  those  that  recurve,  were  associated  with 
large  scale  mid-latitude  tropospheric  flow  anomalies. 

E.  TELECONNECTION  MODELING 

The  anomalies  identified  by  Horel  and  Wallace  (1981),  Wallace  and  Gutzler 
(1981),  and  others  have  many  similarities  with  the  response  to  tropical  forcing  found  in 
teleconnection  modeling  studies  (e.g,  Matsuno  1966,  Egger  1977,  Opsteegh  and  van  den 
Dool  1980,  Hoskins  and  Karoly  1981,  Webster  1981).  In  particular,  Hoskins  and  Karoly 
(1981)  suggested  that  the  arcing  pattern  of  Rossby  wave  energy  propagation  from 
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the  tropics  into  the  mid-latitudes  found  in  many  modeling  studies  might  explain  arcing 
patterns  such  as  the  PNA  pattern. 

Simmons  et  al.  (1983)  found  that  a  variety  of  localized  forcings  in  the  northern 
hemisphere  tropics  and  subtropics  tended  to  produce  telecoimection  patterns  similar  to 
the  empirical  patterns  foimd  by  Wallace  and  Gutzler  (1981).  Simmons  et  al.  (1983) 
attributed  the  model  teleconnection  patterns  to  the  most  rapidly  growing  mode  associated 
with  the  barotropic  instability  of  a  zonally  variable  basic  state.  In  effect,  the  model 
response  to  the  localized  forcing  grew  via  the  extraction  of  kinetic  energy,  particularly  in 
the  vicinity  of  strong  jets.  Branstator  (1983)  and  Lau  and  Lim  (1985)  suggested  that 
westerly  jets  may  act  as  a  barrier  to  and  guide  for  propagating  upper  tropospheric  waves. 

Short  term  teleconnection  studies  using  tropical  cyclones  as  the  tropical  forcing 
were  conducted  by  Woll  (1993),  Springer  (1994),  and  Jakus  (1995).  These  three  studies 
used  an  operational  global  data  assimilation  system  with  tropical  cyclone  bogusing  and  a 
numerieal  weather  prediction  model.  They  foimd  that  tropical  cyclones  may  produce  a 
strong  response  over  the  North  Pacific  -  North  American  (NPNA)  region  which  persists 
well  after  the  tropical  cyclone  dissipated.  These  studies  also  demonstrated  that: 

•  Tropical  cyclones  can  induce  teleconnections, 

•  Teleconnections  are  strongly  affected  by  wave  guiding  and  amplification 
associated  with  the  east  Asian  -  North  Pacific  jet. 
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•  The  absolute  vorticity  flux  associated  with  the  divergent  outflow  from  the 
tropical  cyclone  crossing  the  jet  may  initiate  the  wave  response  to  the 
storm. 

F.  DESIGN  OF  THIS  STUDY 

This  study  builds  upon  the  previous  work  of  Woll  (1993),  Springer  (1994),  and 
Jakus  (1995)  to  examine  the  response  of  a  global  atmosphere  to  individual  tropical 
cyclones  with  varying  tracks  and  intensities.  This  study  differs  from  these  three  previous 
studies  in  that  it  uses  a  simplified  general  circulation  model.  The  usefulness  of  such  a 
model  for  this  kind  of  teleconnection  study  was  determined  in  part  by  comparison  to  the 
results  of  Woll  (1993),  Springer  (1994),  and  Jakus  (1995). 

Our  intent  was  to  develop  a  model  that  is  inexpensive  to  run  and  easy  to  modify,  so  that 
numerous  case  studies  could  be  conducted. 

The  major  hypotheses  tested  in  this  study  were  based  on  the  results  and 
speculations  produced  in  the  three  case  studies  conducted  by  Woll  (1993),  Springer 
(1994),  and  Jakus  (1995).  These  hypotheses  are  given  below. 

•  Tropical  cyclones  can  produce  strong  teleconnections. 

•  The  development  of  teleconnections  is  strongly  influenced  by  wave 
guiding  and  amplification  effects  associated  with  the  mid-latitude  westerly 
jet. 


7 


•  The  absolute  vorticity  flux  associated  with  the  divergent  outflow  from  the 
tropical  cyclone  crossing  the  jet  may  initiate  the  wave  response  to  the 
storm. 

•  Barotropic  instability  (BTI)  processes  play  a  major  role  in  the 
amplification  process. 

•  The  basic  teleconnection  process  varies  with  season  and  the  forcing 
location  due  to  temporal  and  spatial  variations  in  the  ambient  flow, 
especially  the  westerly  jets. 

•  The  results  of  the  simplified  general  circulation  model  and  the  results  of 
the  NWP  model  are  qualitatively  similar  because:  (1)  the  ambient  upper 
tropospheric  flow  plays  a  dominant  role  in  the  development  of  the 
teleconnection  response;  and  (2)  the  ambient  flow  is  similar  in  the  two 
models. 

The  modeling  and  analysis  methods  used  in  this  study  are  described  in  Chapter  II. 
The  model  results  are  documented  in  Chapter  III.  Conclusions  and  recommendations  for 
further  study  are  presented  in  Chapter  IV. 
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II.  METHODS 


A.  MODEL  DESCRIPTION 

The  atmospheric  model  used  in  this  study  is  a  nonlinear  spectral  shallow  water 
equation  model  (GSWM)  on  a  global  spherical  domain.  The  model  has  basic  state 
forcings  that  maintain  a  realistic  basic  state  circulation.  A  perturbation  forcing  is  used  in 
the  model  to  represent  a  tropical  cyclone  for  which  the  amplitude  and  location  may  vary 
both  temporally  and  spatially. 

B.  MODEL  EQUATIONS 

The  derivation  of  the  basic  set  of  shallow  water  equations  and  the  spectral 
formulations  is  described  by  Bourke  (1972)  and  Haltiner  and  Williams  (1980).  The 
derivation  for  the  equations  used  in  this  study  is  similar  but  includes  dissipation,  basic 
state  forcing,  and  perturbation  forcing  terms.  The  basic  governing  equations  of  the  model 
are  the  equation  of  motion  and  the  continuity  equation: 

(1) 


— i=-$  v-r+a  -A  +e 


(2) 


9 


V  is  the  total  horizontal  wind  vector,  O,  =  ^  with  $  e.=  gh^  being  a  constant  global 

mean  geopotential  with  a  fixed  equivalent  depth  h^  and  $  being  the  deviation  from  e- 
The  Qb  forcing  terms  maintain  zonal  and  meridional  variations  in  the  basic  state  during 
the  model  run.  Qp  represents  the  perturbation  forcing.  The  A  terms  are  the  model's 
dissipation  terms. 

An  expression  for  vorticity  is  foimd  by  taking  the  curl  of  (1).  Divergence  is  found 
by  taking  the  divergence  of  (2).  Applying  the  Helmholtz  theorem,  the  model  atmosphere 
circulation  can  be  represented  in  terms  of  stream  fimction,  T,  and  the  velocity 
potential,x,  with  the  expressions: 

C=V2T  (3) 

(4) 

It  is  more  convenient  to  replace  the  equation  of  motion  with  vorticity  and  divergence 
equations  in  spectral  models  because  their  relation  is  simplified  when  spherical  fimctions 
are  used  as  basis  fimctions.  An  additional  benefit  of  using  T  and  x  is  that  the  equations 
contain  only  scalar  terms  (Haltiner  and  Williams  1980). 

Equations  (3)  and  (4)  are  substituted  into  (1)  and  (2).  The  resulting  equations  are 
then  converted  into  spherical  coordinates.  The  dissipation  terms  in  (1)  and  (2)  are 
expanded  into  damping  and  diffusion  terms.  The  damping  relaxes  the  perturbation  state 
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back  to  the  basic  state  on  a  time  scale  of  1 5  days.  The  diffusion  terms  transfer  vorticity, 
divergence,  or  potential  energy  from  small  to  large  length  scales  with  diffusivity 
coefficients:  Kv,  K^,  and  Kg .  Each  diffusivity  coefficient  has  a  value  equal  to  10'^  mVs. 

To  solve  the  spherical  coordinate  forms  of  the  basic  shallow  water  equations,  they 
were  spectrally  formulated  with  a  rhomboidal  25  (R25)  truncation  that  corresponds  to  a 
spherical  Gaussian  grid  resolution  of  2.8“  latitude  by  4.7“  longitude.  The  resulting  model 
equations  are: 

dX  dip 

-2Q  (sintp  + 

(C/T) 

-2Q  (sin(p  V^Y-C/^/a) -V2[(C/^'+F^0/2  cos^  (p  +  O] 

d^/dt  =  -1  /  (ocos^(p)  — (t/,^)+cos<p  — (TO) 

[5a  5(p 


=  -1  /  (acos^  (p)  ^(T  V2’P)+cos  (p-^ 
dt  [  dX  dip 
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The  horizontal  wind  components  are  represented  by: 


u  =  U^l  cos  (p  = 


-(cos<p) 


d(p 


dx 


/(flcosq)) 


(8) 


V  =  V^/  cos  (p  = 


9Y  ,  ,dx 

—  +(COS(p)^ 
dA  d(p 


/  (acoscp) 


(9) 


A.  and  (p  longitude  and  latitude  coordinates.  Q  is  Earth's  angular  velocity  and  a  is  Earth's 
radius.  These  equations  are  solved  via  semi-implicit  differencing  and  the  transform 
method  (Haltiner  and  Williams  1980).  Implicit  and  explicit  leapfrog  steps  were  used 
with  0.5  hour  time  steps  for  all  model  runs.  The  model  fields,  Y,  x,  and  were 
initialized  with  the  basic  state  fields  used  for  each  particular  model  run. 

The  model  equations  describe  the  horizontal  circulation  and  height  fields  of  the 
model  atmosphere  for  a  single  vertical  mode.  The  equivalent  depth,  h^,  determines  the 
vertical  mode.  For  this  study,  an  equivalent  depth  of  h*  =  2000  m  was  used  for  all  model 
runs.  This  depth  was  chosen  to  represent  the  observed  low  frequency  global  mean 
tropospheric  structure,  accounting  for  both  the  first  internal  baroclinic  mode  of  the 
tropical  troposphere  (Gill  1980)  and  the  external  equivalent  barotropic  mode  of  the 
extratropics  (Sawyer  1970,  Wallace  1987). 
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1.  Basic  state  Forcing 

Qb  terms  were  calculated  for  each  basic  state  and  then  held  constant  through  the 
model  runs.  These  terms  are  the  negatives  of  the  residuals  of  the  steady  state  forms  of  the 
expanded  spherical  forms  of  the  vorticity  and  divergence  equations  in  the  absence  of 
perturbation  forcing.  The  basic  state  forcing  terms  are  used  to  maintain  a  zonally  and 
meridionally  varying  ambient  atmospheric.  These  terms  give  a  rough  approximation  of 
orographic  effects,  land-sea  contrast  effects,  radiation,  and  other  processes  not  explicitly 
contained  in  the  model. 

The  basic  state  fields  used  in  this  study  were  the  European  Center  for  Mediiun- 
range  Forecasting  (ECMWF)  analyses  of  global  long  term  monthly  mean  values  of  200 
mb  stream  function,  velocity  potential,  and  geopotential  height.  For  most  of  the  runs  in 
this  study,  the  January  and  July  fields  were  used.  These  two  months  showed  the  best 
contrast  of  summer  and  winter  tropospheric  conditions. 

Figure  2  shows  the  January  basic  state  fields.  A  prominent  feature  is  the  strong 
east  Asian  -  North  Pacific  jet  with  a  jet  maximum  southeast  of  Japan  (Figure  2a).  Other 
features  of  the  January  basic  state  include  a  region  of  diffluence  over  the  western  coast  of 
North  America  and  a  weak  North  American  subtropical  jet  over  the  southeastern  U.S. 
coast  (Figure  2a).  The  southern  hemisphere  200  mb  basic  state  is  much  more  zonal  in 
nature. 
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In  the  July  basic  state  (Figure  3),  the  east  Asian  -  North  Pacific  jet  and  the  North 
American  jet  are  weaker  and  have  shifted  poleward,  and  there  is  a  wide  band  of  tropical 
easterlies  encircling  the  globe  in  the  northern  tropics.  The  southern  hemisphere  exhibits  a 
very  zonal  structure  with  a  strong  jet  over  Australia  and  the  southwest  Pacific. 

Additional  features  of  the  January  and  July  basic  states  are  discussed  in  the  Analysis 
Methods  section  of  this  chapter. 

2.  Perturbation  Forcing 

The  equation  for  the  spatial  variation  of  the  perturbation  forcing,  Qp ,  is: 

0^(X,(p)=0^^exp[-(a/25)(cos'*(cos(A-A,p)7C/M^xcos((p-(Pp)7i:/iV^))^]  (10) 

The  term  s  is  the  perturbation  forcing  scale  length  in  meters.  Xq  and  <po  are  the  central 
longitude  and  latitude  of  the  perturbation  forcing.  Mr  and  Nr  are  the  total  number  of 
model  gridpoints  in  the  longitudinal  and  latitudinal  directions,  respectively. 

Qpo  is  the  amplitude  factor  for  the  perturbation  forcing.  This  term  defines  the 
central  and  maximum  magnitude  of  the  perturbation  forcing.  This  value,  with  imits  of  m’ 

^  s'^ ,  may  be  converted  via  the  hypsometric  equation  to  an  equivalent  heating  rate.  In  this 
study,  the  maximum  perturbation  forcing  used  was  equivalent  to  a  heating  rate  of 
20K/day  or  less. 

The  perturbation  forcing  is  intended  to  represent  a  strong  and  relatively  long- 
lasting  tropical  cyclone.  Thus,  the  parameters  in  (5)  were  selected  to  produce  a  forcing 
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for  which  the  region  with  equivalent  heating  of  IK/day  or  more  has  a  spatial  scale  of 
about  2000  km  or  less.  The  forcing  is  allowed  to  vary  in  its  amplitude  and  location,  with 
the  amplitude  variations  being  linear  functions  of  time.  The  amplitude  and  spatial  scale 
of  the  forcing  are  representative  of  observed  and  theoretical  tropical  cyclones  of  typhoon 
strength  (Elsberry  1985).  The  location  variations  simulate  the  tracks  of  recixrving  and 
straight  running  tropical  cyclones,  with  the  recurving  forcing  extending  into  the  tropics 
and  the  straight  running  forcing  staying  within  28°  of  the  equator.  For  most  of  the  runs  in 
this  study,  the  duration  of  the  forcing  is  the  first  11.5  days  of  the  run.  Some  basic 
features  of  the  perturbation  forcing  are  shown  in  Figure  4. 

C.  MODEL  OUTPUT  FIELDS 

The  model  output  fields  used  in  this  study  have  a  spatial  resolution  of  2.8  degrees 
latitude  by  4.7  degrees  longitude  and  a  temporal  resolution  of  12  hours.  The  follo\ving 
fields  were  examined.  Because  the  basic  state  fields  are  based  on  analyses  at  200  mb,  the 
model  output  fields,  and  the  fields  derived  from  them  (e.g.,  wave  activity  flux)  are 
representations  of  the  atmosphere  at  200  mb. 

1.  Geopotential  Height 

The  geopotential  field  is  especially  effective  in  showing  the  model's  response  to 
tropical  forcing  and  in  identifying  teleconnection  mechanisms. 
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2.  Stream  Function 


The  model's  stream  function  was  used  to  calculate  the  rotational  component  of  the 
horizontal  wind.  This  field  and  the  model's  velocity  potential  field  were  used  to  calculate 
the  model's  total  wind  field. 

3.  Velocity  Potential 

The  model’s  velocity  potential  was  used  to  calculate  the  divergent  component  of 
the  horizontal  wind.  The  divergent  winds  are  capable  of  advecting  absolute  vorticity  and 
are  key  component  of  Rossby  wave  source  calculations  (Sardeshmukh  and  Hoskins 
1988). 

D.  MODEL  RUNS 

The  standard  model  runs  simulate  a  24-day  period  with  constant  basic  state 
forcing  throughout  the  run  and  perturbation  forcing  only  during  the  first  1 1.5  days.  We 
conducted  several  model  runs,  outlined  in  Table  1,  to  study  the  response  of  the  model 
atmosphere  to  short  term  localized  forcing.  The  primary  model  parameters  that  we  varied 
for  the  different  runs  were  the  perturbation  forcing's  amplitude,  Qp^,  and  location.  We 
used  different  basic  states  for  the  runs  but  focused  on  the  January  and  July  basic  states. 

The  model's  equivalent  depth,  perturbation  forcing  scale  length,  dissipation  terms,  and 
other  parameters  were  not  varied  for  this  study.  Only  one  perturbation  forcing  was  used 
for  each  run.  We  conducted  three  categories  of  model  runs  which  are  described  below. 
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1.  Basic  Runs 


The  basic  runs  were  the  simplest  runs  used  in  this  study.  These  runs  were 
conducted  for  both  January  and  July  basic  states  and  with  the  perturbation  forcing  being 
located  in  the  tropical  western  North  Pacific.  These  runs  were  conducted  to  establish  the 
baseline  responses  to  tropical  forcing. 

In  the  first  set  of  these  runs,  the  perturbation  forcing  was  held  at  a  constant 
amplitude  and  location  for  the  first  1 1 .5  days  of  the  model  run  then  turned  off  for  the  last 
12.5  days  of  the  model  run.  The  fixed  amplitude  and  location  were  the  averages  of  those 
used  in  the  recurving  and  straight  running  cases  described  earlier. 

In  the  second  set  of  basic  runs,  the  location  of  the  perturbation  forcing  was  fixed 
as  in  the  first  set  of  basic  runs.  But  the  forcing  amplitude  was  linearly  increased  during 
the  first  six  days  and  linearly  decreased  during  the  next  six  days  to  simulate  the 
strengthening  and  weakening  of  a  tropical  cyclone. 

In  the  third  set  of  basic  runs,  the  perturbation  forcing  location  was  varied  while  its 
amplitude  was  held  constant  as  in  the  first  set  of  basic  runs.  The  location  variations 
simulated  recurving  and  straight  running  cases. 

2.  Tropical  Basin  Runs 

The  tropical  basin  runs  reproduced  the  first  set  of  the  basic  runs,  with  fixed 
perturbation  forcing  amplitude  and  location  but  with  the  forcing  located  in  three  tropical 
ocean  basins:  Bay  of  Bengal,  the  western  North  Pacific,  and  the  western  North  Atlantic. 
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These  runs  were  done  to  establish  the  basic  sensitivity  of  the  response  to  the  tropical 
basin  on  which  the  forcing  was  located. 

3.  Tropical  Cyclone  Runs 

The  tropical  cyclone  runs  give  the  most  realistic  simulation  of  tropical  cyclone 
forcing.  In  these  runs,  the  full  variability  of  the  perturbation  forcing  is  used.  Thus  this 
forcing  had  variable  amplitude,  recurving  and  straight  miming  tracks,  and  was  located  in 
all  three  basins.  For  the  nms  using  different  basins,  the  amplitudes  and  tracks  of  the 
perturbation  forcing  were  the  same.  These  runs  were  done  to  test  the  sensitivity  to  a 
combination  of  forcing  changes. 

E.  ANALYSIS  METHODS 

The  model  response  to  the  perturbation  forcing  is  defined  for  this  study  as  the 
difference  between  a  model  run  with  perturbation  forcing  and  a  comparable  run  without 
perturbation  forcing,  the  difference  being  the  former  minus  the  latter.  This  response  was 
analyzed  using  several  diagnostic  tools.  The  basic  response  patterns  were  identified 
using  the  geopotential  height  fields. 

1.  Quasi-Geostrophic  Wave  Activity  Flux  Vectors 

Plumb  (1985)  developed  a  method  of  tracking  queisi-geostrophic  (QG)  wave 
activity  fluxes  based  on  conservation  relations  and  the  phjise  independent  nature  of  QG 
wave  activity.  The  flux  vectors  he  developed  are  parallel  to  the  path  of  wave  energy 
propagation  while  the  vector  magnitude  are  proportional  to  the  strength  of  the 
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propagation.  Divergent  flux  vectors  are  indicative  of  a  source  of  QG  wave  energy,  while 
convergent  flux  vectors  indicate  a  wave  energy  sink.  Plumb  flux  vectors  and  flux 
divergence  were  not  calculated  within  10°  of  the  equator  due  to  the  inherent  weakness  of 
the  QG  approximation  in  the  tropics. 

2.  Barotropic  Instability  (BTI) 

As  discussed  in  Chapter  I,  barotropic  instability  processes  may  be  involved  in  the 
development  of  teleconnection  responses.  In  this  case,  energy  from  the  background  flow 
is  directed  into  the  response  to  produce  an  amplified  response.  Areas  of  potential  BTI  are 
strongly  associated  with  the  flanks  of  the  extratropical  jets.  Areas  of  potential  BTI  were 
calculated  using  the  Rayleigh-Kuo  instability  criterion  (Kuo  1949): 


P 


dy^ 


=0 


(11) 


P  is  the  total  derivative  of  the  Coriolis  parameter.  The  Rayleigh-Kuo  criterion  is  a 
necessary,  but  not  sufficient  condition  for  the  existence  of  barotropic  instability. 
Therefore,  the  Rayleigh-Kuo  criterion  only  indicates  areas  where  BTI  might  be  involved 
in  the  teleconnection  process.  The  model's  zonal  wind  fields,  U200,  and  p  were  used  to 
calculate  regions  where  (1 1)  was  satisfied.  Figures  2e  and  3e  show  the  areas  of  potential 
BTI  for  the  January  and  July  basic  states.  Note  in  particular  the  large  regions  along  the 
jets  (see  Figures  2a  and  3a)  where  (1 1)  is  satisfied. 
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3.  Stationary  Wave  Number 

To  estimate  the  impacts  of  the  basic  state  flow  on  the  meridional  propagation  of 
low  frequency  QG  wave  energy,  the  stationary  wave  number  (Karoly  1978)  was 
calculated  for  the  January  and  July  basic  states.  This  quantity  gives  an  estimate  of  the 
largest  zonal  wave  number  for  which  meridional  propagation  is  possible,  and  is  inversely 
proportional  to  the  zonal  wind: 

(12) 

u 

Here,  K,  is  the  stationary  wave  number,  K,i  is  the  internal  Rossby  radius  of  deformation  = 
(fr/ghj)'^,  and  u  is  the  zonal  flow.  This  relation  shows  that  strong  zonal  flows  limit 
meridional  propagation,  with  the  limitation  being  greater  for  waves  with  smaller  zonal 
scales.  Figures  2f  and  3f  show  for  the  January  and  July  basic  states.  In  these 

figures,  each  contour  indicates  the  limit  of  meridional  propagation  for  waves  with  zonal 
wave  numbers  less  than  or  equal  to  that  designated  by  the  contour.  Note  that  the  K^  =  5 
contour  tends  to  coincide  with  the  equatorward  flanks  of  the  extratropical  jets  (see 
Figures  2a  and  3a).  This  indicates  that  the  jets  may  be  a  significant  wave  guide  for  waves 
with  zonal  wave  numbers  ^  5. 
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4.  Rossby  Wave  Source  Calculation 

Sardeshmukh  and  Hoskins  (1988)  suggested  that  the  advection  of  absolute 
vorticity  by  the  horizontal  component  of  the  divergent  wind  could  be  a  significant  process 
for  generating  Rossby  waves.  Previous  teleconnection  modeling  studies  ignored  the 
divergent  portion  of  the  horizontal  wind  and  focused  solely  on  the  rotational  component 
of  the  horizontal  wind,  due  to  its  greater  magnitude  on  the  large  scales  associated  with 
teleconnections.  Sardeshmukh  and  Hoskins  (1988)  concluded  that  a  tropical  forcing  may 
not  directly  generate  an  extratropical  Rossby  wave  trains,  but  rather  do  so  by  producing 
strong  divergent  winds  that  cross  regions  of  large  absolute  vorticity  gradients  along  the 
jets.  This  in  turn  would  lead  to  a  large  wave  generation  term  in  the  absolute  vorticity 
equation. 

The  isobaric  coordinate  form  of  the  absolute  vorticity  equation  contains  the 
Rossby  wave  source  term,  S: 


Where  (C+f)  is  the  absolute  vorticity  term. 


(13) 


21 


The  Rossby  source  term,  S,  is  defined  as: 

-5=-V-[K^(C+/)] 


(14) 


Expansion  of  (14)  yields: 


s=- 


^  0«„  dv.  ^ 


[  dx  dy}\ 


(15) 


The  first  three  terms  on  the  right-hand  side  of  (15)  represent  advection  of  absolute 
vorticity  by  the  divergent  wind  (ADVDIV)  .  The  other  source  term  in  (15)  is  the  vortex 
stretching  or  divergence  term. 

The  terms  in  (15)  were  calculated  for  various  times  in  the  model  runs  using  finite 
difference  approximations  of  the  derivatives.  The  differences  between  these  source  terms 
and  those  in  the  initial  fields  represent  the  Rossby  wave  source  response  to  the 
perturbation  forcing. 
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III.  RESULTS 


A.  INTRODUCTION  TO  MODEL  RUNS 

To  test  the  model's  response  to  short  term  forcing,  a  series  of  model  runs  was 
developed  and  executed.  These  model  runs  began  with  basic  forcing  scenarios  and 
progressed  to  model  runs  that  simulated  the  more  complex  forcing  behavior  associated 
with  tropical  cyclones  (Table  1). 

Before  conducting  the  model  experiments,  we  performed  a  series  of  model  tests  to 
determine  the  basic  characteristics  and  range  of  responses  to  different  basic  state 
conditions  and  perturbation  forcing  values.  These  tests  showed  that  the  results  from  runs 
with  the  January  (July)  basic  state  conditions  were  very  representative  of  the  results  from 
runs  with  other  northern  winter  (summer)  basic  state  conditions.  In  addition,  these  tests 
showed  that  the  results  obtained  when  using  a  January  basic  state  or  a  July  basic  state 
were  markedly  different  from  each  other.  Thus,  we  selected  as  our  primary  basic  state 
conditions  those  from  January  and  July  in  order  to  highlight  the  effects  of  seasonal 
variations  in  the  basic  state  using  a  minimum  of  runs. 

Our  tests  of  the  model's  sensitivity  to  the  forcing  amplitude  showed  qualitatively 
similar  responses  for  all  reasonable  values  of  the  amplitude.  Thus,  we  used  forcing 
amplitudes  that  were  20  K/day,  with  the  amplitude  being  constant  for  some  runs  and 
varying  with  time  for  other  runs  (see  Table  1). 
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Recurving  and  straight  running  forcing  paths  were  developed  for  the  variable 
forcing  location  runs.  These  paths  were  based  upon  the  characteristic  tracks  of  actual 
tropical  cyclones.  The  recurving  path  (Figure  4a)  allowed  the  perturbation  forcing  to 
penetrate  the  mid-latitudes,  while  the  straight  running  path  remained  entirely  in  the 
tropics  below  28°  N  (Figure  4b). 

In  the  following  sections  of  this  chapter,  we  review  the  major  results  from  the 
runs  conducted  for  this  study  and  summarized  in  Table  1.  Our  focus  in  this  review  is  on 
the  northern  hemisphere  responses  to  the  tropical  perturbation  forcing.  This  is  because: 
(1)  our  forcing  experiments  all  used  northern  hemisphere  perturbation  forcing  and  gave 
responses  that  were  relatively  large  in  the  northern  hemisphere;  and  (2)  we  wanted  to 
compare  our  results  to  those  from  the  studies  of  Woll  (1993),  Springer  (1994),  and  Jakus 
(1995)  which  focused  on  responses  in  the  North  Pacific  -  North  American  region. 

B.  WESTERN  NORTH  PACIFIC  RUN  DESCRIPTIONS 

1.  January  Fixed  Forcing  Runs 

Figure  5  shows  the  fixed  forcing  and  height  response  for  the  most  basic  January 
western  North  Pacific  (run  PJl,  Table  1).  The  response  after  12  hours  (Figure  5a)  was  a 
positive  height  response  over  southern  Japan  north-northwest  of  the  perturbation  forcing 
location.  A  wave  pattern  developed  rapidly  across  the  North  Pacific  to  the  east  of  this 
early  response.  This  wave  pattern  had  a  very  slow  westward  drift  during  the  first  few 
days. 
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However,  by  day  four  of  the  model  run,  the  initial  positive  response  became  fixed  over 
the  west  coast  of  China  (Figure  5b). 

This  fixed  pattern  continued  to  amplify  eastward  from  the  east  Asian  region  along 
approximately  35®  N.  The  pattern  encircled  the  northern  hemisphere  by  day  eight  with  a 
clear  zonal  wavenumber  four  component  (Figure  5c).  During  most  of  the  run,  the 
amplitude  in  the  initial  response  region  decreased  as  the  amplitude  increased 
progressively  further  to  the  east.  This  suggests  an  eastward  group  velocity  for  the  wave 
response.  The  response  persisted  well  after  the  perturbation  forcing  ceased  (Figure  5b-g). 
Throughout  most  of  the  nm,  the  high  and  low  height  responses  over  the  northeast  Pacific 
had  a  pronounced  northwest-southeast  tilt  and  roughly  resembled  the  negative  of  the  PNA 
pattern. 

a.  Rossby  Wave  Source  Analyses 

Figure  6  shows  the  early  responses  for  some  of  the  terms  in  the  Rossby 
wave  source  equation,  (1 5).  The  divergent  wind  response  (Figure  6a)  showed  a  distinct 
radial  outflow  from  the  perturbation  forcing  region,  with  strong  flow  crossing  the  jet  over 
east  Asia  (Figure  2a).  The  ADVDIV  component  of  the  Rossby  wave  source  response 
(Figure  6b)  had  a  distinct  negative  feature  just  south  of  Japan  that  coincided  with  the 
position  of  the  initial  positive  height  response  (Figure  5a).  This  ADVDIV  response 
persisted  throughout  the  first  12  days  of  the  run  (not  shown).  The  stretching  component 
of  the  Rossby  wave  source  response  (Figure  6c)  had  a  negative  feature  over  the 
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perturbation  forcing  region  which  persisted  throughout  the  forcing  period  (the  first  1 1.5 
days)(not  shown).  There  were  also  distinct  positive  and  negative  features  over  northern 
Japan  and  northeastern  Asia.  These  features  were  associated  with  weak  height  responses 
(e.g.,  the  positive  height  response  southeast  of  Lake  Baikal,  Figure  5a).  Throughout  the 
first  12  days,  these  stretching  responses  varied  considerably  and  were  associated  with 
relatively  weak  and  transient  height  responses  (not  shown). 

Throughout  the  first  12  days  of  the  model  run,  the  ADVDIV  and  stretching 
responses  had  the  same  magnitude,  but  only  the  negative  ADVDIV  response  was 
associated  with  the  positive  height  response  over  east  Asia.  The  negative  source  response 
is  consistent  with  the  positive  height  response,  since  both  are  associated  with  the 
generation  of  negative  relative  vorticity.  Thus,  it  appears  that  the  ADVDIV  response 
may  have  helped  initiate  the  height  response. 

b.  Wave  Energy  Analyses 

The  QG  wave  activity  flux  vectors  corresponding  to  the  height  responses 
in  Figure  5  are  shown  in  Figure  7.  There  was  a  weak  energy  propagation  away  from  the 
forcing  region  at  12  hours  (Figure  7a).  The  region  of  strongest  propagation  shifted 
eastward  during  the  first  12  days  until  it  became  centered  over  the  northeast  Pacific 
(Figure  7a-d).  This  was  also  a  region  of  strong  northeastward  propagation.  To  the  west 
of  this  region,  the  propagation  was  predominantly  eastwcird  along  the  east  Asian  -  North 
Pacific  jet  (e.g.,  Figure  7c).  Over  North  America,  the  propagation  described  an  arcing 
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pattern,  with  the  energy  directed  southeastward  over  central  Canada  (e.g.,  Figure  7d). 

The  divergence  of  the  flux  vectors  in  Figure  8  indicates  strong  energy  sources  for  the 
wave  response  over  the  subtropical  North  Pacific.  Note,  for  example,  the  divergent 
vectors  near  Hawaii  after  eight  days  (Figure  7c),  on  the  south  flank  of  the  basic  state  jet 
near  the  jet  exit  (see  Figure  2a). 

The  strong  zonal  propagation  along  the  jet  over  the  North  Pacific,  and  the 
strong  meridional  propagation  in  the  midlatitude  inter-jet  regions  over  the  North  Pacific 
and  North  America,  suggest  that  the  jet  acted  as  a  wave  guide  (cf  Lau  and  Lim  1984). 
This  is  consistent  with  the  restricted  meridional  propagation  along  the  jets  indicated  by 
the  basic  state  stationary  wave  numbers  (Figure  2f).  The  strong  wave  energy  source 
along  the  North  Pacific  jet  flank  suggests  that  the  jet  was  an  energy  source  for  the 
response.  The  major  energy  source  region  in  the  subtropical  North  Pacific  is  collocated 
with  the  downstream  end  of  the  region  of  potential  BTI  in  the  basic  state  (Figure  2e). 

This  suggests  that  BTI  processes  were  involved  in  developing  this  wave  energy  source. 
Similar  results  were  found  by  Woll  (1993)  and  Springer  (1994). 

2.  January  Fixed  Location- Variable  Strength  Runs 

The  response  patterns  for  the  January  fixed  location-varied  amplitude  run  (run 
PJl  A,  Table  1)  were  qualitatively  similar  to  those  for  the  fixed  forcing  run  (run  PJl, 
Table  1  and  preceding  section).  The  height  response  after  12  days  (Figure  8)  is  similar  to 
that  in  run  PJl  except  that  it  is  weaker  and  has  a  more  zonal  structure  over  the  western 
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North  Pacific  (cf.  Figure  5d).  The  Rossby  wave  sources  are  also  similar  for  the  two  runs. 
The  energy  propagation  and  energy  sources  (not  shown)  are  also  similar  to  but  more 
zonal  than  those  for  run  PJl  (Figure  7). 

3.  January  Variable  Location-Fixed  Strength  Runs 

Figure  9  shows  the  height  response  to  the  western  North  Pacific  run  with 
recurving,  fixed  amplitude  forcing  (run  PJIL,  Table  1)  after  12  days.  This  response,  and 
the  Rossby  wave  source  and  wave  energy  responses  (not  shown),  were  similar  to  but 
markedly  stronger  than  those  for  runs  PJl  and  PJl  A  (see  table  1  and  the  preceding 
sections  of  this  chapter).  This  was  especially  true  when  the  forcing  approached  the  east 
Asian  jet.  There  was  also  a  significant  energy  source  on  the  poleward  flank  of  the  east 
Asian  -  North  Pacific  jet  (not  shown)  in  the  region  of  potential  BTI  in  the  basic  state 
(Figure  2e).  The  movement  of  the  perturbation  forcing  into  the  extratropics  may  explain 
the  development  of  this  midlatitude  energy  source  which  was  not  so  apparent  in  the  runs 
for  which  the  forcing  location  was  fixed.  These  differences  in  the  runs  suggest  that  the 
same  basic  mechanism  was  operating  in  all  three  cases  and  that  interactions  between  the 
perturbation  forcing  and  the  jet  (e.g.,  Rossby  wave  source  processes)  may  be  significant 
in  producing  the  wave  response. 

4.  January  Recurving  Cyclone  Run 

The  responses  to  the  recurving,  variable  amplitude  forcing  (run  PJIR,  Table  1) 
were  similar  to  that  for  run  PJIL  but  with  weaker  responses  (cf  Figures  9, 10  for  the 
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height  responses).  This  difference  is  consistent  with  a  weaker  average  forcing  in  run 

PJIR,  especially  when  the  forcing  entered  the  extratropics.  The  energy  propagation 
response  showed  strong  wave  guiding  along  the  jet  in  the  western  North  Pacific  during 
days  1-10  (not  shown).  During  days  12-18,  there  was  a  strong  arcing  energy  propagation 
path  out  of  the  northeastern  Pacific  and  across  North  America  at  12  days  (Figure  10). 

The  major  source  region  for  this  energy  was  the  downstream  end  of  the  basic  state  region 
of  potential  BTI  over  the  subtropical  North  Pacific  (Figure  2a). 

5.  January  Straight  Running  Cyclone  Run 

The  responses  to  the  straight  running,  variable  amplitude  tropical  forcing  (nm 

PJIS,  Table  1)  were  very  similar  to  but  weaker  than  those  for  the  recurving  variation  of 
this  nm  (nm  PJIR),  as  shown  by  a  comparison  of  Figures  10  and  11.  These  results  are 
consistent  with  the  previously  discussed  runs,  and  support  the  conclusion  that  the 
processes  by  which  the  response  develops  are  relatively  insensitive  to  the  location  of  the 
forcing  --  at  least  so  long  as  the  divergent  outflow  from  the  forcing  reaches  the  east  Asian 
-  North  Pacific  jet. 

6.  July  Fixed  Forcing  Runs 

The  height  response  to  the  run  using  a  fixed  forcing  in  the  western  Pacific  and  a 
July  basic  state  (run  PJ2,  Table  1)  is  shown  in  Figure  12.  At  12  hours  into  the  run  there 
was  a  positive  height  response  over  the  perturbation  forcing  (Figure  12a).  By  four  days, 
this  positive  response  had  shifted  about  1000  km  to  the  west  of  the  forcing,  and  a  weak 
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zonal  wave  train  had  developed  along  approximately  40®N  (Figure  12b).  This  wave  train 
had  a  clear  zonal  wave  number  three  component  and  a  westward  phase  propagation  of 
about  13  m/s.  During  the  next  eight  days,  this  basic  response  pattern  remained,  with  the 
midlatitude  response  growing  in  both  the  northern  and  southern  hemispheres  (Figure 
12c,d).  After  the  perturbation  forcing  was  turned  off  after  1 1 .5  days,  the  positive 
response  near  the  forcing  quickly  disappeared,  but  the  wave  train  at  40°N  persisted  and 
maintained  its  strength  throughout  most  of  the  run.  This  was  especially  true  over  east 
Asia,  where  the  waves  amplified  as  they  propagated  westward  through  the  region. 
Throughout  the  run,  the  northern  hemisphere  response  was  as  strong  or  stronger  than  that 
in  the  southern  hemisphere. 

a.  Rossby  Wave  Source  Analyses 

Figure  13  shows  the  early  responses  for  some  of  the  terms  in  the  Rossby 
wave  source  equation,  (15).  The  divergent  wind  response  (Figure  13a)  showed  a  distinct 
radial  outflow  from  the  perturbation  forcing  region,  with  strong  flow  crossing  the  jet  over 
east  Asia  (Figure  3a).  The  ADVDIV  component  of  the  Rossby  wave  source  response 
(Figure  1 3b)  had  a  weak  negative  feature  over  the  Sea  of  Japan  that  coincided  with  the 
position  of  a  weak  positive  height  response  (Figure  12a).  The  stretching  component  of 
the  Rossby  wave  source  response  (Figure  13  c)  had  a  strong  negative  feature  over  the 
perturbation  forcing  region  which  coincided  with  the  positive  height  response  there 
(Figure  12a).  The  stretching  response  also  had  distinct  positive  and  negative  features 
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over  northeastern  Asia,  but  these  did  not  coincide  with  any  significant  height  responses. 
The  collocation  over  the  Sea  of  Japan  of  the  weak  negative  ADVDIV  response  and  the 
weak  positive  height  response  at  12  hours  suggests  that  the  ADVDIV  process  may  have 
made  a  small  contribution  to  the  initiation  of  the  response  in  this  region.  However,  the 
propagation  of  positive  and  negative  height  responses  through  this  region  suggests  that 
the  ADVDIV  process  was  not  so  important  in  determining  the  response  as  was  the  case 
for  the  eorresponding  January  run,  PJl.  This  difference  was  probably  due  to  the  weaker 
east  Asian  jet  in  the  July  basic  state. 

b.  Wave  Energy  Analyses 

The  QG  wave  activity  flux  vectors  corresponding  to  the  height  responses 
in  Figure  12  are  shown  in  Figure  14.  During  most  of  the  first  12  days,  there  was  energy 
propagation  to  the  northeast,  east,  and/or  southeast  away  from  the  forcing  region  (Figure 
14a-d),  indieating  that  the  perturbation  forcing  was  a  significant  source  of  wave  energy. 
This  is  consistent  with  the  strong  stretching  response  seen  in  this  area  during  the  first  12 
days  (e.g..  Figure  13c).  During  the  remaining  days,  the  major  energy  source  was  over  the 
east  Asian  -  northwest  Pacific  region  at  about  40-50°N.  Figure  3e  indicates  that  this 
region  was  potentially  barotropieally  unstable,  suggesting  that  the  basic  state  flow  may 
have  contributed  energy  to  the  response  in  this  region.  From  this  region,  the  energy 
propagated  eastward,  opposite  to  the  phase  propagation  indicated  by  the  height  responses 
(Figure  12).  The  zonal  energy  propagation  during  days  12-24  across  east  Asia  and 
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the  northwest  Pacific  at  about  45°N  is  consistent  with  the  zonal  wave  guide  indicated  by 
the  stationary  wave  numbers  in  this  area  (Figure  3f). 

The  results  from  this  run,  PJ2,  contrast  sharply  with  those  from  the 
corresponding  run  which  used  a  January  basic,  PJl  (see  Chapter  III,  section  Bl).  These 
differences  are  due  to  the  differences  in  the  basic  states,  since  those  are  the  only 
differences  between  the  runs  (compare  Figures  2  and  3).  The  weaker  (stronger)  northern 
(southern)  hemisphere  responses  in  PJ2  indicate  that  the  strength  of  the  ambient  flow  is  a 
significant  factor  in  allowing  a  large  response  to  develop. 

7.  July  Fixed  Location-Variable  Strength  Runs 

The  July  fixed  location,  variable  amplitude  run  (PJ2A,  Table  1)  response  (not 
shown)  was  qualitatively  similar  to  that  for  run  PJ2,  but  much  weaker  in  the  northern 
hemisphere.  This  suggests  that  a  significant  northern  hemisphere  response  in  the 
presence  of  a  weak  basic  state  flow  requires  a  relatively  sustained  tropical  forcing  or, 
possibly,  a  forcing  that  travels  closer  to  the  extratropics. 

8.  July  Variable  Location-Fixed  Strength  Runs 

The  response  after  12  and  16  days  of  the  July  run  in  which  the  forcing  recurved 
into  the  extratropics  with  a  fixed  amplitude  (run  PJ2L,  Table  1)  is  shown  in  Figure  15. 
The  positive  response  over  the  Sea  of  Japan  in  Figure  1 5a  was  a  feature  that  traveled 
northwestward  with  the  forcing  during  the  first  12  days.  This  feature  dissipated  within  a 
few  days  of  the  forcing  being  turned  off  after  1 1.5  days.  By  day  16,  the  northern 
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hemisphere  response  resembled  that  seen  in  run  PJ2  (Figure  12e).  The  phase 
propagation,  Rossby  wave  source,  and  wave  energy  analyses  (not  shown)  also  indicate 
that  these  two  runs  were  similar.  The  major  differences  between  these  runs  seem  to  be 
the  result  of  the  northwestward  shift  of  the  positive  height  response  over  the  forcing  that 
occurred  in  PJ2L  but  not  in  PJ2. 

9.  July  Recurving  Cyclone  Run 

The  response  of  the  July  run  with  recurving  forcing  with  variable  amplitude  (run 
PJ2R,  Table  1)  was  very  similar  to  that  for  runs  PJ2  and  PJ2L.  The  height  response  after 
12  and  16  days  (Figure  16)  is  similar  to  the  response  at  about  12-20  days  seen  in  PJ2  and 
PJ2L  (Figures  12  and  15).  In  particular,  the  response  in  the  three  runs  had  similar  spatial 
and  temporal  scales,  westward  phase  propagation,  eastward  group  propagation,  and 
amplification  in  areas  of  potential  BTI  in  the  east  Asian  -  North  Pacific  region. 

10.  July  Straight  Running  Cyclone  Run 

The  response  for  run  PJ2S  (not  shown)  was  very  similar  to  that  for  run  PJ2R.  The 
Rossby  wave  source  and  wave  energy  analyses  (not  shown)  that  this  similarity  was  due  to 
the  ability  of  the  forcing  in  the  two  runs  to  generate  a  response  in  the  area  of  potential 
BTI  over  east  Asia  at  40-50°N  (Figure  3e). 

C.  BAY  OF  BENGAL  RUN  DESCRIPTIONS 

The  runs  for  which  the  tropical  perturbation  forcing  was  located  in  the  Bay  of 
Bengal  gave  results  that  were  similar  in  many  ways  to  the  corresponding  runs  with  the 
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forcing  in  the  tropical  western  North  Pacific.  In  particular,  the  Rossby  wave  source  and 
wave  energy  analyses  showed  that  similar  mechanisms  were  operating  for  both  the  Bay 
of  Bengal  and  Pacific  runs.  Because  of  these  similarities,  only  the  key  results  from  some 
of  the  Bay  of  Bengal  runs  are  described  here. 

1.  January  Fixed  Forcing  Runs. 

Figure  17  shows  the  fixed  forcing  and  height  response  for  the  most  basic  January 
Bay  of  Bengal  run  (run  BJl,  Table  1).  The  response  after  12  hours  (Figure  17a)  was  a 
positive  height  response  over  the  Tibetan  Plateau  north  of  the  perturbation  forcing 
location.  A  wave  pattern  developed  rapidly  across  Asia  and  the  North  Pacific  to  the  east 
of  this  early  response.  This  wave  pattern  had  a  very  slow  westward  drift  during  the  first 
few  days.  However,  by  day  four  of  the  model  run,  the  initial  positive  response  became 
fixed  over  the  southwest  China  (Figure  17b). 

This  fixed  pattern  continued  to  amplify  eastward  firom  the  central  Asian  region 
along  approximately  35°  N.  The  pattern  encircled  the  northern  hemisphere  by  day  eight 
with  a  clear  zonal  wavenumber  four  component  (Figure  17c).  During  most  of  the  middle 
and  later  periods  of  the  run,  the  amplitude  in  the  initial  response  region  decreased  as  the 
amplitude  increased  progressively  further  to  the  east.  This  suggests  an  eastward  group 
velocity  for  the  wave  response.  The  response  persisted  well  after  the  perturbation  forcing 
ceased  (Figure  17b-g). 
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Throughout  most  of  the  run,  the  high  and  low  height  responses  over  the  northeast  Pacific 
had  a  pronounced  northwest-southeast  tilt  and  roughly  resembled  the  PNA  pattern. 

These  BJl  responses  were  very  similar  to  those  for  PJl,  the  corresponding 
western  Pacific  run,  except  that  the  response  over  Asia  was  better  developed  in  BJl  and 
the  overall  responses  were  out  of  phase  over  most  of  the  northern  hemisphere.  These 
results  suggest  that  tropical  disturbances  in  the  Indian  Ocean  may  be  just  as  effective  at 
producing  anomalies  in  the  North  Pacific  -  North  American  region  as  those  in  the  tropical 
Pacific. 

a.  Rossby  Wave  Source  Analyses 

A  negative  Rossby  wave  source  response  (not  shown)  developed  early  in 
the  run  at  the  location  of  the  positive  response  over  the  Tibetan  plateau.  At  this  location, 
the  dominant  term  for  this  response  was  the  ADVDIV  term  (not  shown).  Thus,  it  appears 
that  the  ADVDIV  response  may  have  helped  initiate  the  height  response,  as  was  the  case 
for  the  corresponding  western  Pacific  run  (PJl). 

b.  Wave  Energy  Analyses 

The  QG  wave  activity  flux  vectors  corresponding  to  the  height  responses 
in  Figure  17  are  shown  in  Figure  18.  There  was  a  weak  energy  propagation  away  from 
the  forcing  region  at  12  hours  (Figure  18a).  The  region  of  strongest  propagation  shifted 
eastward  during  the  first  12  days  until  it  became  centered  over  the  northeast  Pacific 
(Figure  18a-d).  This  was  also  a  region  of  strong  northeastward  propagation.  To  the  west 
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of  this  region,  the  propagation  was  predominantly  eastward  along  the  east  Asian  -  North 
Pacific  jet  (e.g.,  Figure  1 8c).  In  the  later  part  of  the  run,  the  propagation  over  North 
America  described  a  rough  arcing  pattern,  with  the  energy  directed  southeastward  over 
eastern  Canada  (e.g..  Figure  18e).  The  divergence  of  the  flux  vectors  in  Figure  18 
indicates  strong  energy  sources  for  the  wave  response  over  eastern  Asia  and  the 
subtropical  North  Pacific.  Note,  for  example,  the  divergent  vectors  near  northeastern 
China  and  Hawaii  after  eight  days  (Figure  18c). 

The  strong  zonal  propagation  along  the  jet  over  east  Asia  and  the  North 
Pacific,  and  the  strong  meridional  propagation  in  the  midlatitude  inter-jet  regions  over  the 
North  Pacific  and  North  America,  suggest  that  the  jet  acted  as  a  wave  guide  for  the 
response.  This  is  consistent  with  the  restricted  meridional  propagation  along  the  jets 
indicated  by  the  basic  state  stationary  wave  numbers  (Figure  2f).  The  strong  wave  energy 
sources  along  the  east  Asian  -  North  Pacific  jet  flanks  suggests  that  the  jet  was  an  energy 
source  for  the  response.  Both  of  the  major  energy  sources  were  collocated  with  regions 
of  potential  BTI  in  the  basic  state  (Figure  2e).  This  suggests  that  BTI  processes  were 
involved  in  developing  these  wave  energy  sources.  Similar  results  were  found  for  the 
corresponding  January  run,  P J 1 . 

2.  January  Recurving  Cyclone  Run 

Figure  19  shows  the  height  response  to  the  recurving,  variable  amplitude  forcing 
after  12  days  (run  BJIR,  Table  1).  The  similarity  between  this  response  and  that  for  run 
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BJl  (Figure  5d)  occurred  at  other  times  and  in  the  other  model  fields  analyzed  (not 
shown).  In  particular,  the  energy  propagation  response  showed  strong  wave  guiding 
along  the  jet  in  the  east  Asian  -  North  Pacific  during  days  1-10  (not  shown).  During  days 
12-18,  there  was  a  strong  arcing  energy  propagation  path  out  of  the  northeastern  Pacific 
and  across  North  America  at  12  days  (Figure  20).  The  major  source  region  for  this 
energy  was  the  downstream  end  of  the  basic  state  region  of  potential  BTI  over  the 
subtropical  North  Pacific  (Figure  2e). 

3.  July  Recurving  Cyclone  Run 

The  response  of  the  July  run  with  recurving  forcing  with  variable  amplitude  (run 
BJ2R,  Table  1)  was  qualitatively  similar  to  that  for  the  corresponding  western  Pacific 
runs  PJ2R  (compare,  for  example.  Figures  21  and  16).  The  BJ2R  and  PJ2R  responses 
had  similar  spatial  and  temporal  scales,  westward  phase  propagation,  eastward  group 
propagation,  and  amplification  in  areas  of  potential  BTI  in  the  east  Asian  -  North  Pacific 
region.  The  main  differences  between  the  two  runs  were  due  to  longitudinal  shifts  in  the 
response  locations. 

D.  NORTH  ATLANTIC  RUN  DESCRIPTIONS 

The  runs  for  which  the  tropical  perturbation  forcing  was  located  in  the  North 
Atlantic  gave  results  that  had  some  clear  similarities  and  differences  when  compared  with 
the  results  from  the  corresponding  tropical  western  North  Pacific  or  Bay  of  Bengal  runs 
(see  Table  1). 
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The  key  results  from  some  of  the  North  Atlantic  runs  are  described  here  to  highlight  these 
similarities  and  differences. 

1.  January  Fixed  Forcing  Runs 

Figure  22  shows  the  fixed  forcing  and  height  response  for  the  most  basic  January 
North  Atlantic  run  (run  All,  Table  1).  The  response  after  12  hours  (Figure  22a)  was  a 
positive  height  response  approximately  collocated  with  the  perturbation  forcing  northeast 
of  Colombia  (Figure  22a).  A  wave  pattern  developed  rapidly  from  North  America 
eastward  across  the  North  Atlantic  and  into  Asia  and  the  North  Pacific  to  the  east  of  this 
early  response  (Figure  22a-d).  This  wave  pattern  had  a  very  slow  westward  drift  during 
the  first  few  days.  However,  by  day  eight  of  the  model  run,  the  positive  and  negative 
centers  of  the  response  were  approximately  fixed  in  place  (Figure  22c). 

This  fixed  pattern  amplified  eastward  from  the  North  Atlantic  region  along 
approximately  30-50“N.  This  suggests  an  eastward  group  velocity  for  the  wave  response. 
The  pattern  encircled  the  northern  hemisphere  by  day  eight  with  a  clear  zonal 
wavenumber  four  component  (Figure  22c).  Throughout  the  run,  the  strongest  responses 
were  in  the  North  Atlantic  -  European  region.  The  response  persisted  well  after  the 
perturbation  forcing  ceased  (Figure  22e-g).  Throughout  most  of  the  run,  the  high  and  low 
height  responses  over  the  North  Atlantic  -  European  region  had  a  pronounced  northwest- 
southeast  tilt. 
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This  tilt  was  also  seen  during  the  middle  part  of  the  run  in  the  North  Pacific  -  North 
American  region  where  a  pattern  roughly  similar  to  the  PNA  pattern  developed  (Figure 
22d-e). 

The  overall  AJl  height  responses  over  Asia  and  the  North  Pacific  were  roughly 
similar  to  those  for  PJl  and  BJl  (Figures  5  and  17).  However,  over  the  North  Atlantic  - 
European  region,  the  AJl  responses  were  stronger,  more  persistent,  more  tilted.  These 
results  suggest  that  tropical  disturbances  in  the  North  Atlantic  region  may  affect  the 
circulation  throughout  the  northern  hemisphere,  but  have  their  greatest  impact  in  the 
North  Atlantic  -  European  area. 

a.  Rossby  Wave  Source  Analyses 

After  12  hours  of  the  AJl  run,  a  weak  negative  ADVDIV  response  (Figure 
23b)  had  developed  in  association  with  the  divergent  wind  response  (Figure  23a)  across 
the  North  Atlantic  jet  (Figure  2a).  However,  the  development  of  the  early  height 
response  was  mainly  associated  with  a  strong  negative  stretching  term  in  the  Rossby 
wave  source  response  (Figure  23c).  This  is  in  contrast  to  the  corresponding  Pacific  and 
Bay  of  Bengal  results  (runs  PJl  and  BJl)  in  which  the  ADVDIV  term  was  better 
correlated  with  the  height  anomalies.  This  difference  is  probably  due  to  the  North 
Atlantic  jet  being  considerably  weaker  than  the  Asian  -  North  Pacific  jet. 
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b.  Wave  Energy  Analyses 

The  QG  wave  activity  flux  vectors  corresponding  to  the  height  responses 
in  Figure  23  are  shown  in  Figure  24.  There  was  a  weak  eastward  and  northeastward 
energy  propagation  away  from  the  forcing  region  at  12  hours  which  grew  stronger  during 
the  first  10  days  (Figure  24a-c).  The  region  of  strongest  propagation  shifted  eastward 
during  the  first  16  days  until  it  became  centered  over  the  northeast  Pacific  (Figure  24a-e). 
This  was  also  a  region  of  northeastward  propagation.  To  the  west  of  this  region,  the 
propagation  was  predominantly  eastward  along  the  north  African  -  south  Asian  -  North 
Pacific  jet  (compare  Figures  24d  and  2a).  In  the  later  part  of  the  run,  the  propagation  over 
the  North  Pacific  -  North  American  region  described  a  weak  arcing  pattern,  with  the 
energy  directed  southeastward  over  the  central  U.S.  (e.g..  Figure  24f).  The  divergence  of 
the  flux  vectors  in  Figure  24  indicates  strong  energy  sources  for  the  wave  response  over 
the  western  North  Atlantic  -  Iberian  -  Mediterranean  region  and  the  subtropical  North 
Pacific.  Note,  for  example,  the  divergent  vectors  near  the  western  Mediterranean  and 
west  of  Hawaii  after  twelve  days  (Figme  24d). 

The  strong  propagation  along  the  jets  over  the  North  Atlantic,  north 
Africa,  southern  Asia,  and  the  North  Pacific,  plus  the  strong  meridional  propagation  in 
the  midlatitude  inter-jet  regions  over  the  North  Pacific,  North  America,  and  North 
Atlantic,  suggest  that  the  jets  acted  as  a  wave  guide  for  the  response.  As  an  example, 
note  the  northeastward  propagation  at  8  and  12  days  between  the  North  Atlantic  jet  and 
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the  north  African  jet  (Figure  24c-d).  The  basic  state  stationary  wave  number  analysis 
(Figure  2f)  indicates  that  this  is  a  region  of  relatively  unrestricted  meridional  propagation, 
or  a  propagation  duct,  between  the  barriers  posed  by  the  two  jets. 

In  addition,  the  strong  wave  energy  sources  along  the  flanks  of  the  North 
Atlantic  jet  and  the  north  African  -  south  Asian  -  North  Pacific  jet  suggests  that  these  jets 
were  energy  sources  for  the  response.  The  major  energy  sources  along  the  jets  were 
collocated  with  regions  of  potential  BTI  in  the  basic  state  (Figure  2e).  This  suggests  that 
BTI  processes  were  involved  in  developing  these  wave  energy  sources.  Similar  results 
were  found  for  the  corresponding  January  run  in  the  Pacific,  PJl,  and  Bay  of  Bengal, 

BJl. 

2.  January  Recurving  Cyclone  Run 

Figure  25  shows  the  height  response  to  the  recurving,  variable  amplitude  forcing 
after  12  days  (run  AJIR,  Table  1).  The  similarity  between  this  response  and  that  for  run 
AJl  (Figure  22d)  occurred  at  other  times  and  in  the  other  model  fields  analyzed  (not 
shown).  In  particular,  there  was  a  strong  energy  source  in  the  North  Atlantic  -  Iberian  - 
Mediterranean  region  (e.g..  Figure  26)  associated  with  the  basic  state  region  of  potential 
BTI  between  the  North  Atlantic  and  north  African  jets  (Figure  2e).  During  the  middle 
parts  of  the  AJIR  run,  the  responses  were  stronger  than  for  AJl,  probably  because  during 
this  period  the  forcing  was  stronger  and  located  closer  to  the  areas  of  potential  BTI. 
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3.  July  Recurving  Cyclone  Run 

The  response  of  the  July  run  with  recurving  forcing  with  variable  amplitude  (run 
AJ2R,  Table  1)  was  similar  to  that  for  the  other  North  Atlantic  July  runs.  After  4  days, 
the  height  response  for  run  AJ2R  (Figure  27b)  had  a  positive  response  over  and  to  the 
northwest  of  the  forcing  and  a  negative  response  to  the  northeast.  This  pattern  persisted 
until  the  forcing  was  turned  off  after  11.5  days  (Figure  27b-d).  During  days  10-20,  an 
arcing  wave  train  occurred  over  the  North  American  -  North  Atlantic  -  European  area 
which  drifted  westward  after  the  forcing  was  turned  off  (Figure  27d-f).  The  response 
over  most  of  Asia  was  noticeably  weaker  and  less  persistent  than  in  the  Pacific  and  Bay 
of  Bengal  July  runs. 

a.  Rossby  Wave  Source  Analyses 

After  12  hours  of  the  AJ2R  run,  a  weak  negative  ADVDIV  response  (not 
shown)  had  developed  in  association  with  the  divergent  wind  response  (not  shown)  across 
the  subtropical  North  Atlantic  jet  (Figure  3a).  However,  the  development  of  the  early 
height  response  was  mainly  associated  with  a  strong  negative  stretching  term  in  the 
Rossby  wave  source  response  (not  shown).  This  was  similar  to  the  corresponding 
January  runs,  AJl  and  AJIR. 

b.  Wave  Energy  Analyses 

The  QG  wave  activity  flux  vectors  corresponding  to  the  height  responses 
at  4,  8,  12,  and  16  days  (Figure  27)  are  shown  in  Figure  28.  There  was  a  weak  eastward 
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and  northeastward  energy  propagation  away  from  the  forcing  region  at  4  days.  At  8  days, 
when  the  forcing  was  near  its  maximum  amplitude  and  located  over  Florida,  there  was 
very  strong  and  divergent  southeastward  and  eastward  propagation  from  the  forcing 
region.  This  energy  source  was  located  at  the  upstream  end  of  an  area  of  potential  BTI 
between  the  North  American  jet  and  the  subtropical  North  Atlantic  jet  (Figure  3a).  By  12 
days,  this  energy  source  was  gone.  During  days  14-20,  there  was  moderate  to  weak 
propagation  on  an  arcing  path  across  North  America.  The  upstream  end  of  this  path  was 
collocated  with  the  downstream  end  of  an  area  of  potential  BTI  on  the  southern  flank  of 
the  North  Pacific  jet  (Figure  3e).  At  this  time  near  12  hours  which  grew  stronger  during 
the  first  10  days  (Figure  24a-c). 

E.  SUMMARY 

The  model  results  showed  a  number  of  patterns  which  are  summarized  in  the 
following  sections. 

1.  Height  Responses 

All  the  runs  showed  clear  teleconnection  responses  to  the  tropical  perturbation 
forcing.  These  responses  were  most  clear  in  the  height  responses,  in  which  distinct  wave 
trains  tended  to  develop.  The  runs  with  relatively  simple  representations  of  a  tropical 
perturbation  forcing  gave  responses  that  were  qualitatively  similar  to  those  foimd  in  the 
corresponding  nms  with  more  realistically  varying  forcing.  This  indicates  that  this  model 
is  relatively  insensitive  to  many  of  the  specifics  of  the  perturbation  forcing.  The  height 
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response  patterns  in  the  July  and  January  runs  were  distinctly  different.  In  particular,  the 
northern  hemisphere  response  in  January  tended  to  be  stronger  and  more  stationary  than 
in  July.  For  the  July  runs,  the  response  had  a  clear  westward  phase  propagation  and  the 
southern  hemisphere  response  was  generally  as  strong  as  the  northern  hemisphere 
response,  even  when  the  forcing  was  located  well  north  of  the  equator.  This  tendency  for 
a  strong  winter  hemisphere  response  has  been  found  in  other  studies  (e.g.,  Opsteegh  and 
Van  den  Dool,  1980). 

There  was  a  tendency  for  the  height  responses  to  start  with  a  positive  response 
somewhere  to  the  north  of  the  forcing  (most  of  the  January  runs)  or  approximately  over 
the  forcing  (the  July  and  North  Atlantic  runs).  The  response  then  developed  by  eastward 
energy  propagation  away  from  the  initial  positive  response.  This  eastward  propagation  of 
the  response  tended  to  be  most  rapid  for  the  North  Atlantic  runs  and  least  rapid  for  the 
Bay  of  Bengal  runs.  The  poleward  propagation  was  also  most  rapid  for  the  North 
Atlantic  runs  and  least  rapid  for  forcing  in  the  Bay  of  Bengal.  This  is  consistent  with  the 
southwest  to  northeast  tilt  of  the  jets  and  the  propagation  ducts  in  the  northeast  Atlantic 
region  (Figures  2a,f  and  3a,f).  The  responses  downwind  of  the  forcing  longitude  tended 
to  be  zonally  shorter  than  those  on  upwind  side.  This  is  consistent  with  the  dispersion 
relationship  for  Rossby  waves  (Holton,  1992).  The  response  patterns  tended  to  shift 
longitudinally  as  the  forcing  was  shifted  between  the  tropical  basins. 
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The  January  runs  generally  had  northeastward  energy  propagation  over  the 
northeast  Pacific  and  the  northeast  Atlantic.  This  is  consistent  with  the  low  wind  speed 
regions  between  the  jets  and  the  propagation  ducts  over  the  northeast  Pacific  -  North 
American  region  and  over  the  northeast  Atlantic  region  (Figures  2a, f  and  3a,f).  The 
January  runs  also  had  response  patterns  that  persisted  well  after  the  forcing  was  turned 
off,  with  the  shorter  waves  decaying  more  rapidly  after  the  forcing  was  turned  off 

The  July  runs  had  response  patterns  that  generally  drifted  to  the  west,  with  the 
response  to  North  Atlantic  forcing  drifting  least  rapidly.  The  responses  tended  to  amplify 
as  they  passed  through  the  region  of  potential  BTI  over  east  Asian  near  northern  Japan. 

The  runs  with  forcing  in  the  Bay  of  Bengal  or  western  North  Pacific  had  similar 
responses,  although  they  tended  to  be  out  of  phase  with  each  other.  The  runs  with  North 
Atlantic  forcing  had  distinctly  different  response  patterns.  However,  for  all  of  these  runs, 
there  tended  to  be  a  clear  response  in  the  northeast  Pacific  -  North  American  region,  with 
the  positive  and  negative  height  response  centers  describing  a  poleward  and  eastward 
arcing  pattern. 

2.  Teleconnection  Dynamics 

fl.  Rossby  Wave  Source  Analyses 

The  initial  positive  height  responses  were  associated  with  a  Rossby  wave 
source  process.  For  the  January  runs  with  forcing  in  the  Bay  of  Bengal  or  in  the  western 
North  Pacific,  this  process  depended  on  the  advection  of  absolute  vorticity  as  the 
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forcing's  divergent  outflow  crossed  a  nearby  subtropical  jet.  For  the  January  runs  with 
North  Atlantic  forcing,  this  process  depended  on  vortex  stretching  within  the  forcing.  For 
the  July  runs,  the  process  depended  mainly  on  vortex  stretching.  These  differences  in  the 
Rossby  wave  initiation  were  due  to  longitudinal  and  seasonal  variations  in  the  position 
and  strength  of  the  jet  poleward  of  the  forcing,  with  the  jet  being:  (1)  stronger  and  closer 
to  the  forcing  in  the  January  Bay  of  Bengal  and  western  North  Pacific  cases;  and  (2) 
weaker  and  further  from  the  equator  in  the  North  Atlantic  and  July  cases. 

b.  Wave  Energy  Analyses 

The  wave  activity  fluxes  confirmed  the  predominantly  eastward  energy 
propagation  indicated  by  the  development  of  height  responses.  These  fluxes  also  showed 
regions  of  strong  poleward  and  equatorward  energy  propagation  in  the  weak  flow  areas 
between  the  jets.  The  basic  state  stationary  wave  number  calculations  indicated  that  these 
areas  were  also  potential  propagation  ducts. 

The  divergence  of  the  wave  activity  fluxes  showed  that  the  major  sources  of 
energy  for  the  response  were:  (1)  the  pertiubation  forcing  during  the  initial  development 
of  the  response;  (2)  areas  of  potential  barotropic  instability  in  the  basic  state.  The 
imstable  area  on  the  southern  flank  of  the  North  Pacific  jet  was  collocated  with  relatively 
large  energy  sources  in  many  of  the  responses,  especially  for  the  January  runs.  This  may 
explain  why  the  northeast  Pacific  -  North  American  region  tended  to  have  large  and 
persistent  responses. 
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IV.  CONCLUSIONS 


A.  SUMMARY 

We  used  a  global  shallow  water  model  to  investigate  the  dynamics  of 
teleconnections  initiated  by  relatively  small  scale,  short  term  tropical  perturbations  (e.g., 
tropical  cyclones).  Our  emphasis  was  on  the  northern  hemisphere  response  to 
perturbations  in  the  northern  tropics.  Our  study  is  a  companion  to  those  by  Woll  (1993), 
Springer  (1994),  and  Jakus  (1995)  who  used  a  global  data  assimilation  system  and 
numerical  weather  prediction  (NWP)  model  to  investigate  the  response  to  individual 
tropical  cyclones  in  the  western  Pacific.  The  model  for  our  study  was  much  easier  to 
modify  and  more  economical  to  run  than  a  NWP  model.  This  allowed  us  to  conduct 
many  runs  and  investigate  the  importance  of  many  variables  (e.g.,  variations  in  the  path, 
amplitude,  and  longitude  of  the  perturbation  forcing;  the  seasonal  variations  of  the  basic 
state). 

The  major  hypotheses  tested  in  this  study,  and  in  those  by  Woll  (1993),  Springer 
(1994),  and  Jakus  (1995),  are  given  below. 

•  Tropical  cyclones  can  produce  strong  teleconnections. 

•  The  development  of  teleconnections  is  strongly  influenced  by  wave 
guiding  and  amplification  effects  associated  with  the  mid-latitude  westerly 
jet. 


47 


•  The  absolute  vorticity  flux  associated  with  the  divergent  outflow  from  the 
tropical  cyclone  crossing  the  jet  may  initiate  the  wave  response  to  the 
storm. 

•  •  Barotropic  instability  (BTI)  processes  play  a  major  role  in  the 

amplification  process. 

•  The  basic  teleconnection  process  varies  with  season  and  the  forcing 
location  due  to  temporal  and  spatial  variations  in  the  ambient  flow, 
especially  the  westerly  jets. 

•  The  results  of  the  simplified  general  circulation  model  and  the  results  of 
the  NWP  model  are  qualitatively  similar  because:  (1)  the  ambient  upper 
tropospheric  flow  plays  a  dominant  role  in  the  development  of  the 
teleconnection  response;  and  (2)  the  ambient  flow  is  similar  in  the  two 
models. 

The  results  of  this  study  confirm  all  of  these  hypotheses.  In  particular,  they 
indicate  that  the  simplified  global  model  used  in  this  study  may  yield  results 
that  are  at  least  qualitatively  consistent  with  those  from  a  much  more  sophisticated  NWP 
model. 

B.  IMPLICATIONS  FOR  EXTENDED-RANGE  FORECASTING 

The  upper  tropospheric  height  responses  seen  in  this  and  previous  studies  (e.g., 
Woll  1993,  Springer  1994,  Jakus  1995,  and  Harr  and  Elsberry  1995)  indicate  that 
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relatively  small  scale  and  short  term  tropical  disturbances,  such  as  tropical  cyclones,  may 
have  significant  global  impacts  on  time  scales  of  one  to  three  weeks.  The  mechanism 
results  from  this  study  suggest  that  these  impacts  are  more  likely  at  certain  times  and 
locations.  For  example,  strong  and  persistent  impacts  are  more  likely  when  the  tropical 
disturbance  is  able  to  interact  with  a  strong  jet.  Thus,  when  the  disturbance  travels  well 
away  fi'om  the  equator,  or  when  the  jet  shifts  toward  the  equator,  significant 
teleconnections  are  more  likely.  Thus,  in  the  development  of  extended  range  forecast 
models,  close  attention  should  be  paid  to  accurately  simulating  such  tropical  disturbances 
and  their  potential  interactions  with  the  jets.  In  addition,  special  attention  should  be  paid 
to  locations  and  times  when  such  interactions  are  more  likely  (e.g.,  in  the  tropical  western 
North  Pacific  during  northern  hemisphere  winter  when  recurving  cyclones  may  readily 
approach  the  subtropical  east  Asian-North  Pacific  jet). 

C.  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

We  recommend  additional  studies  with  the  global  shallow  water  model  because  of 
the  similarity  of  its  results  to  those  from  more  sophisticated  models  and  its  convenience 
and  low  cost.  We  would,  however,  recommend  using  a  considerably  higher  resolution 
version  of  the  model  (e.g.,  R  80).  This  would  help  in  producing  tropical  forcings  with 
more  realistic  areal  extents  and  tracks,  and  basic  state  fields  with  more  realistic 
waveguides  and  instability  features.  This  increased  realism  would  help  in  assessing  the 
relative  importance  of  variations  in  the  forcing  and  interactions  of  the  forcing  with  the 
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basic  state.  In  addition  to  more  studies  with  the  simple  global  model,  we  also 
recommend  additional  NWP  modeling  studies,  combined  with  diagnostic  studies  of 
observed  teleconnections. 

Some  of  the  key  cases  we  recommend  for  these  studies  are: 

1.  Studies  in  which  the  basic  state  is  more  realistic  (e.g.,  a  basic  state  calculated  as 
the  mean  of  the  analyzed  conditions  during  a  particular  two  week  period;  a  time 
varying  basic  state) . 

2.  Studies  using  more  realistic  tropical  disturbances  (e.g.,  multiple  simultaneous 
typhoons  in  the  same  and/or  different  basins  and  possible  interference  effects; 
multiple  sequential  disturbances  and  their  cumulative  impacts) . 

3.  NWP  studies  using  straight  running  tropical  cyclones  and  tropical  cyclones  in 
the  Indian  and  Atlantic  basins. 

4.  Studies  in  which  the  analytic  tools  applied  to  the  model  results  are  applied  to 
observed  anomalies  (e.g.,  the  anomalies  used  by  Harr  and  Elsberry  1995). 

5.  Studies  in  which  the  energy  tendencies  and  conversions  involved  in  the 
development  of  short  term  teleconnections  are  examined  further. 

6.  Studies  investigating  the  impacts  of  short  term  teleconnection  periods  on  the 
midlatitude  atmosphere  and  ocean  (e.g.,  on  the  circulation  of  the  northeast  Pacific 
ocean). 
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Such  studies  should  increase  our  understanding  of  short  term  teleconnections  and 
the  improvement  of  extended-range  foreeasting  teehniques.  The  dynamical  links  between 
short  term  and  long  term  teleeonnections  indicates  that  the  above  studies  may  also  help 
improve  our  understanding  seasonal  and  longer  period  teleconnections  and  the  dynamies 
of  global  elimate  ehange. 
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or  July);  and  (c)  the  location  and  amplitude  ot  the  perturbation  forcing  (fixed  or  variable).  The  run  names  are  indicated 
by  the  three-  or  four-character  names  shown  in  the  cells  of  the  table.  The  first  character  indicates  the  basin  (e.g.,  P  for 
Pacific);  the  second  two  indicate  the  basic  state  month  (e.g.,  J1  for  January);  and  the  fourth  indicates  the  variable  factor 
in  the  perturbation  forcing  (e.g.,  A  for  amplitude,  R  for  recurving  location). 


APPENDIX  A  -  FIGURES 
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Figure  1 .  The  Pacific-North  American  (PNA)  teleconnection  pattern.  The  hypothesized 
global  pattern  of  tropospheric  height  anomalies  (solid  lines)  during  a  northern  hemisphere 
winter  which  coincides  with  an  episode  of  anomalously  warm  sea  surface  temperatures  in 
the  central  and  eastern  equatorial  Pacific.  The  bold  arrows  in  dark  type  represent 
strengthened  subtropical  jets.  The  light  arrows  depict  streamline  distortions  caused  by  the 
height  anomaly  pattern  of  troughing  over  the  central  Pacific  and  ridging  over  western 
Canada.  The  shaded  area  represents  enhanced  cloudiness  and  rainfall.  (From  Horel  and 
Wallace  1981) 


b.  Geopotential  height 


Figure  2.  200  mb  January  basic  state  fields.  Solid  (dashed)  contours  indicate  positive 
(negative)  values,  (a)  Zonal  wind  component,  contour  interval  is  5  m/s.  (b)  Geopotential 
heights  after  subtracting  a  global  mean  value  of  1 1,784  gpm,  contour  interval  is  100  gpm. 
(c)  Streamftinction,  contour  interval  is  1x10’  mVs.  (d)  Velocity  potential,  contour  interval 
is  1x10*  mVs. 


Figure  2.  (Continued)  (e)  Indications  of  potential  for  barotropic  instability.  Zero  contour 
indicates  regions  where  P-Uyy=0.  (f)  Stationary  wavenumber  with  each  contour  indicating 
the  limit  of  meridional  propagation  for  low  frequency  quasi-geostrophic  waves  with  zonal 
wavenumber  less  than  or  equal  to  the  number  on  the  contour,  for  wave  numbers  ^  8. 
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Figure  3.  200  mb  July  basic  state  fields.  Solid  (dashed)  contours  indicate  positive 
(negative)  values,  (a)  Zonal  wind  component,  contour  interval  is  5  m/s.  (b)  Geopotential 
heights  after  subtracting  a  global  mean  value  of  1 1,784  gpm,  contour  interval  is  100  gpm. 
(c)  Streamfunction,  contour  interval  is  1x10’  mVs.(d)  Velocity  potential,  contour  interval 
is  1x10®  mVs. 
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Figure  3.  (Continued)  (e)  Indications  of  potential  for  barotropic  instability.  Zero  contour 
indicates  regions  where  p-Uyy=0.  (f)  Stationary  wavenumber  with  each  contour  indicating 
the  limit  of  meridional  propagation  for  low  frequency  quasi-geostrophic  waves  with  zonal 
wavenumber  less  than  or  equal  to  the  number  on  the  contour,  for  wave  munbers  ^  8. 
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Figure  4.  The  model's  tropical  perturbation  forcing  for  two  western  North  Pacific 
examples.  The  small  oval  surrounds  the  region  in  which  the  forcing  is  equivalent  to 
heating  of  s  IK/day.  The  oval  is  located  at  the  initial  position,  while  the  line  leading 
northwestward  from  the  center  of  the  oval  shows  the  forcing’s  successive  positions  in:  (a) 
the  recurving  case  and  (b)  the  straight  running  case.  The  dots  along  the  path  indicate  the 
forcing's  central  positions  at  48  hour  intervals.  The  most  poleward  dot  shows  the  position 
at  1 1 .5  days,  after  which  the  forcing  was  turned  off.  For  the  Bay  of  Bengal  and  North 
Atlantic  runs,  the  forcing  was  simply  shifted  to  the  west  or  east,  with  no  other  changes. 
For  a  summary  of  all  the  runs  and  their  different  forcing  conditions,  see  Table  1. 
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Figure  5.  Height  response  for  run  PJl  which  used:  (1)  a  January  basic  state;  and  (2)  a 
fixed  location,  fixed  amplitude  tropical  perturbation  forcing  in  the  western  North  Pacific. 
For  each  panel,  the  elapsed  model  time  is  indicated.  The  small  oval  contour  in  the  tropical 
western  North  Pacific  surrounds  the  region  in  which  the  forcing  was  equivalent  to  heating 
^  1  K/day.  Positive  (negative)  responses  are  indicated  by  the  solid  (dashed)  contours. 
Contour  interval  is  2  gpm.  Zero  contour  omitted. 
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Figure  5.  (Continued) 


f.  20  days 


Figure  5.  (Continued) 


g.  24  days 


Figure  5.  (Continued) 


Figure  6.  200  mb  Rossby  wave  source  response  terms  for  run  PJl  after  12  hours  of  the 
model  run.  (a)  Divergent  wind  response,  with  vector  scale,  in  ms"',  as  shown  below,  (b) 
Advection  by  divergent  wind  (ADVDIV)  response,  (c)  vortex  stretching  response. 
Positive  (negative)  responses  are  indicated  by  the  solid  (dashed)  contours.  Contour 
interval  is  0.5x10  "  s  i  Zero  contour  omitted. 
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Figure  7.  Quasi-geostrophic  wave  activity  flux  response  for  run  PJl .  For  each  panel,  the 
elapsed  model  time  is  indicated.  Tropical  western  North  Pacific  perturbation  forcing 
of  ^  1  K/day  is  indicated  by  the  small  oval  region.  Vector  scale,  in  mVs-  as  shown. 
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Figure  7.  (Continued) 
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Figure  7.  (Continued) 
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12  days 


Figure  8.  Height  response  for  run  PJl  A  which  used:  (1)  a  January  basic  state;  and  (2)  a 
fixed  location,  variable  amplitude  tropical  perturbation  forcing  in  the  western  North 
Pacific.  The  elapsed  model  time  is  indicated.  Positive  (negative)  responses  are  indicated 
by  the  solid  (dashed)  contours.  Contour  interval  is  2  gpm.  Zero  contour  omitted. 
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12  days 


Figure  9.  Height  response  for  run  PJIL  which  used:  (1)  a  January  basic  state;  (2)  a 
recurving  path;  and  (3)  fixed  amplitude  tropical  perturbation  forcing  in  the  western  North 
Pacific.  The  elapsed  model  time  is  indicated.  Positive  (negative)  responses  are  indicated 
by  the  solid  (dashed)  contours.  Contour  interval  is  2  gpm.  Zero  contour  omitted. 
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12  days 


Figure  10.  Height  response  for  run  PJIR  which  used:  (1)  a  January  basic  state;  (2)  a 
recurving  path;  and  (3)  variable  amplitude  tropical  perturbation  forcing  in  the  western 
North  Pacific.  The  elapsed  model  time  is  indicated.  Positive  (negative)  responses  are 
indicated  by  the  solid  (dashed)  contours.  Contour  interval  is  2  gpm.  Zero  contour 
omitted. 
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Figure  11.  Height  response  for  run  PJIS  which  used:  (1)  a  January  basic  state;  (2)  a 
straight  running  path;  and  (3)  variable  amplitude  tropical  perturbation  forcing  in  the 
western  North  Pacific.  The  elapsed  model  time  is  indicated.  Positive  (negative)  responses 
are  indicated  by  the  solid  (dashed)  contours.  Contour  interval  is  2  gpm.  Zero  contour 
omitted. 
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Figure  12.  200  mb  height  response  for  run  PJ2  which  used:  (1)  a  July  basic  state;  and  (2) 
a  fixed  location,  fixed  amplitude  tropical  perturbation  forcing  in  the  western  North 
Pacific.  For  each  panel,  the  elapsed  model  time  is  indicated.  The  small  oval  contour  in  the 
tropical  western  North  Pacific  surrounds  the  region  in  which  the  forcing  was  equivalent 
to  heating  ^  1  K/day.  Positive  (negative)  responses  are  indicated  by  the  solid  (dashed) 
contours.  Contour  interval  is  2  gpm.  Zero  contour  omitted. 
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Figure  12.  (Continued) 
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Figure  12.  (Continued) 
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Figure  12.  (Continued) 
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Figure  13.  200  mb  Rossby  wave  source  response  terms  for  run  PJ2  after  12  hours  of  the 
model  run.  (a)  Divergent  wind  response,  with  vector  scale,  in  ms  ',  as  shown  below,  (b) 
Advection  by  divergent  wind  (ADVDIV)  response,  (c)  vortex  stretching  response. 
Positive  (negative)  responses  are  indicated  by  the  solid  (dashed)  contours.  Contour 
interval  is  0.5x10"  s"'^.  Zero  contour  omitted. 


0  150E*<1^ 
riAx  iHuTTvEmfi 


82 


Figure  14.  Quasi-geostrophic  wave  activity  flux  response  for  run  PJ2.  For  each  panel,  the 
elapsed  model  time  is  indicated.  Tropical  perturbation  forcing  of  ^  1  K/day  is  indicated 
by  the  small  oval  region.  Vector  scale,  in  m^/s^  as  shown. 

Vector  scale,  in  m^/s^  as  shown. 


d.  12  days 


Figure  14.  (Continued) 
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f.  20  days 


Figure  14  (Continued) 
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Figure  15.  Height  response  for  run  PJ2L  which  used:  (1)  a  July  basic  state;  (2)  a 
recurving  path;  and  (3)  fixed  amplitude  tropical  pertubation  forcing  in  the  western  North 
Pacific.  The  elapsed  model  time  is  indicated.  Positive  (negative)  responses  are  indicated 
by  the  solid  (dashed)  contours.  Contour  interval  is  2  gpm.  Zero  contour  is  omitted. 
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Figure  16.  Height  response  for  run  PJ2R  which  used:  (1)  a  July  basic  state;  (2)  a 
recurving  path;  and  (3)  varible  amplitude  tropical  pertubation  forcing  in  the  western 
North  Pacific.  The  elapsed  model  time  is  indicated.  Positive  (negative)  responses  are 
indicated  by  the  solid  (dashed)  contours.  Contour  interval  is  2  gpm.  Zero  contour  is 
omitted. 


Figure  17.  Height  response  for  run  BJl  which  used:  (1)  a  January  basic  state;  (2)  a  fixed 
location;  and  (3)  fixed  amplitude  tropical  pertubation  forcing  in  the  Bay  of  Bengal.  The 
elapsed  model  time  is  indicated.  Positive  (negative)  responses  are  indicated  by  the  solid 
(dashed)  contours.  Contour  interval  is  2  gpm.  Zero  contour  is  omitted. 
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d.  12  days 


Figure  17.  (Continued) 
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f.  20  days 


Figure  17.  (Continued) 
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Figure  17.  (Continued) 
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Figure  1 8.  Quasi-geostrophic  wave  activity  flux  response  for  run  BJl .  For  each  panel, 
the  elapsed  model  time  is  indicated.  Tropical  Bay  of  Bengal  perturbation  forcing 
of  i  1  K/day  is  indicated  by  the  small  oval  region.  Vector  scale,  in  mVs^  as  shown. 


Figure  18.  (Continued) 
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f.  20  days 


Figure  18.  (Continued) 


Figure  19.  Height  response  for  run  BJIR  which  used:  (1)  a  January  basic  state;  (2)  a 
recurving  path;  and  (3)  variable  amplitude  tropical  pertubation  forcing  in  the  Bay  of 
Bengal.  The  elapsed  model  time  is  indicated.  Positive  (negative)  responses  are  indicated 
by  the  solid  (dashed)  contours.  Contour  interval  is  2  gpm.  Zero  contour  is  omitted. 


Figure  20.  Quasi-geostrophic  wave  activity  flux  response  for  run  BJIR.  For  each  panel, 
the  elapsed  model  time  is  indicated.  Tropical  perturbation  forcing  of  ^  1  K/day  is 
indicated  by  the  small  oval  region.  Vector  scale,  in  mVs^  as  shown. 
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Figure  21.  Height  response  for  run  BJ2R  which  used:  (1)  a  July  basic  state;  (2)  a 
recurving  path;  and  (3)  variable  amplitude  tropical  pertubation  forcing  in  the  Bay  of 
Bengal.  The  elapsed  model  time  is  indicated.  Positive  (negative)  responses  are  indicated 
by  the  solid  (dashed)  contours.  Contour  interval  is  2  gpm.  Zero  contoiur  is  omitted. 


Figure  22.  Height  response  for  run  AJl  which  used:  (1)  a  Januaryy  basic  state;  (2)  a  fixed 
location;  and  (3)  fixed  amplitude  tropical  pertubation  forcing  in  the  North  Atlantic.  The 
elapsed  model  time  is  indicated.  Positive  (negative)  responses  are  indicated  by  the  solid 
(dashed)  contours.  Contour  interval  is  2  gpm.  Zero  contour  is  omitted. 
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Figure  22  (Continued) 
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f.  20  days 


Figure  22.  (Continued) 
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Figure  22.  (Continued) 
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Figure  23.  200  mb  Rossby  wave  source  response  terms  for  run  AJl  after  12  hours  of  the 
model  run.  (a)  Divergent  wind  response,  with  vector  scale,  in  ms'',  as  shown  below,  (b) 
Advection  by  divergent  wind  (ADVDIV)  response,  (c)  vortex  stretching  response. 
Positive  (negative)  responses  are  indicated  by  the  solid  (dashed)  contours.  Contour 
interval  is  0.5x10"  s'^.  Zero  contour  omitted. 
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Figure  23  (Continued) 
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Figure  24.  Quasi-geostrophic  wave  activity  flux  response  for  run  AJl .  For  each  panel, 
the  elapsed  model  time  is  indicated.  Tropical  western  North  Pacific  perturbation  forcing 
of  ^  1  K/day  is  indicated  by  the  small  oval  region.  Vector  scale,  in  mVs^  as  shown. 


0.100Et06 
HAXIMUn  VKTOR 


107 


Figure  24.  (Continued) 
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Figure  24.  (Continued) 


a.  12  days 


Figure  25.  Height  response  for  run  AJIR  which  used:  (1)  a  July  basic  state;  (2)  a 
recurving  path;  and  (3)  variable  amplitude  tropical  pertubation  forcing  in  the  North 
Atlantic.  The  elapsed  model  time  is  indicated.  Positive  (negative)  responses  are  indicated 
by  the  solid  (dashed)  contours.  Contour  interval  is  2  gpm.  Zero  contour  is  omitted. 
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Figure  26.  Quasi-geostrophic  wave  activity  flux  response  for  run  AJIR.  For  each  panel, 
the  elapsed  model  time  is  indicated.  Tropical  perturbation  forcing  of  ^  1  K/day  is 
indicated  by  the  small  oval  region.  Vector  scale,  in  mVs^  as  shown. 
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Figure  27.  Height  response  for  run  AJ2R  which  used:  (1)  a  July  basic  state;  (2)  a 
recurving  path;  and  (3)  variable  amplitude  tropical  pertubation  forcing  in  the  North 
Atlantic.  The  elapsed  model  time  is  indicated.  Positive  (negative)  responses  are  indicated 
by  the  solid  (dashed)  contours.  Contour  interval  is  2  gpm.  Zero  contour  is  omitted. 
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Figure  27.  (Continued) 
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Figure  27.  (Continued) 


Figure  27.  (Continued) 
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Figure  28.  Quasi-geostrophic  wave  activity  flux  response  for  run  AJ2R.  For  each  panel, 
the  elapsed  model  time  is  indicated.  Tropical  pertxirbation  forcing  of  ^  1  K/day  is 
indicated  by  the  small  oval  region.  Vector  scale,  in  mVs^  as  shown. 
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d.  16  days 


Figure  28  (Continued) 
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